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Abstract 

We study the effect of fluctuations in opening angle of bent-core molecules on stability of the 
nematic phases. The molecules are built out of two Gay-Berne centers of k = 4, k' = 5, ^ = 1 
and v = 2, each corresponding to one of the molecular arms. Constant-pressure Monte Carlo (MC 
NPT) simulations are carried out for two versions of the model, with fixed- and fluctuating opening 
angle, where fluctuations are taken harmonic about the angle of 140°. The systems studied are first 
cooled down from the isotropic liquid to a highly ordered crystalline smectic phase and then heated 
up all the way back. The liquid crystalline phases found on cooling are uniaxial nematic, biaxial 
nematic and crystalline, hexagonal biaxial smectic. On heating the biaxial nematic is superseded 
by the metastable crystalline phase. Similar system with molecules of shorter arms {k = 3) was 
shown to exhibit a wide range of reduced density, where the uniaxial nematic was absolutely stable 
[Liq. Cryst. 29, 483 (2002)]. For molecules with rigid arms of k = 4 and opening angles of 125°, 
130°, 135° and 140° a further broadening of the density range of the uniaxial nematic phase is 
observed. However, the helical superstructures reported for k = 3 do not show up for k = 4. 

PACS numbers: 61.30.-v, 64.70.Md, 05.10.Ln 



*e-mail address wj503@york.ac.uk 
'''e-mail address |Tech.longa@uj.edu.pI 



I. INTRODUCTION 



The most common phases found for bent-core hquid crystals are smectics [l| . It has there- 
fore been quite surprising to discover that these systems can also stabilize the thermotropic 
uniaxial (Nu)- and biaxial nematic (Nb) phases js, ^, |4j_^, with quest for A''^ occupying 
soft matter research for more than 30 years now (see e.g. To date the factors crucial for 
absolute stability of the Nb phase in bent-core systems are not fully understood. In this 
work we carry out MC NPT simulations for model bent-core (V-shaped) molecules to seek 
for the minimal molecular features necessary to stabilize Nb- An excellent recent review on 
theoretical- and computer simulation attempts to find the Nb phase in model systems is 
found in j3] . A brief summary of the results relevant for what follows is given below. 

The lattice simulations with 'bent-core' dispersion centers occupying lattice sites of a 
cubic lattice show that the Nb phase can be stabilized about the Landau point at, or 
close to, the 'tetrahedral' opening angle of cos' -i(-l/3) = 109.47° fl, y, Q. Additional 
inclusion of the dipole-dipole interactions with dipole moments along the C2 molecular axis 
can transform the Landau point into a Landau line [lo| . 

The presently existing simulations for model bent-core mesogens with centers of mass 
free to move (off-lattice) are less optimistic. The excluded volume effects (packing entropy) 
make the stable nematic phases, especially the biaxial one, rare to appear in these systems. 
For example, an ensemble of bend-core molecules formed by two rigid spherocylinders joined 
at one end stabilize only Njj and isotropic fluid 11, I2I. Replacing spherocylinders with 
Gay-Berne (GB) centers of = 5, /i = 1 and 1/ = 2, also does not bring a stable 



biaxial nematic, as shown by Memmer 



13| and Johnston et al Q, [l^. Even the Nu phase 



disappears for an intermediate opening angle of 170°. Similar conclusions as concerning 
stability of Nu can be drawn from NPT MC simulations of Dewar and Camp [16|] for the 
multi-site bent-core model of rigidly connected LJ spherical sites. Before proceeding further 
we would like to mention about an interesting observation in aforementioned simulations of 
Memmer 13| and Johnston et al [l^. It concerns the appearance of a spontaneously chiral, 
helical state close to the uniaxial nematic-smectic phase transition. Memmer, however, has 
made it clear that the existence of the helical structure might be an artefact of the method 
he used for this structure was not reproduced in every run. 

The Metropolis MC simulations for model bent-core molecules with dipole moments have 
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also been carried out with an attempt to bringing the molecular modeling close to the 
experimental situation. Unfortunately, the presence of a single dipole located on the C2 
molecular axis has led so far to neither Nb nor Nu jl5,[l3, But two point dipoles placed 
at the centers of the lateral sites of a three GB bent-core molecule were shown to stabilize the 



Nu phase in NPT MC simulations 181]. Moving the dipoles towards the terminal positions 



of the molecular arms broadened the uniaxial nematic range. Finally, the molecules with 



terminal flexible chains also do not help to stabilize Nu or Nb [17|, but the simultaneous 
presence of flexible tails and molecular dipoles along the C2 axis can lead to the Nu phase 
17(1 . However, in none of the off-lattice simulations mentioned so far the biaxial nematic 
phase was detected. 

To date, the only off-lattice simulation for bent-core systems that reveals the biaxial 
nematic phase has been that of Pelaez and Wilson 



191 ]. They performed the Molecular 
Dynamics simulation for real molecules with experimentally detected Nb phase usmg 
a full-atomistic potential. The biaxial nematic phase with local ferroelectric domains was 
observed on cooling down an isotropic phase. However, the same atomistic model with 
electrostatic charges switched off gave a stable smectic phase instead of Nb- This seems to 
indicate that a subtle balance between an average bent molecular shape, molecular flexibility 
and transverse dipoles is essential for bent-core mesogenic molecules to stabilize the Nb 
phase, and that the question of the minimal molecular features necessary to obtain Nb is 
still open. 

In this work we present results of Metropolis MC NPT simulations aimed at seeking 
for such minimal molecular features necessary to stabilize Nb for V-shaped molecules. We 
check as whether the elongation of the molecular arms can play any role. We also study 
some consequences of fluctuations in opening angle. We choose the model similar to that 
of Memmer [l^, where bent-core molecules were built out of two Gay-Berne (GB) centers 
representing molecule's arms. An attractive feature of the model 13[ is that it shows the 
stable Nu phase in a relatively wide temperature range. Our model differs from the original 
one by longer molecular arms. In addition we allow for harmonic fluctuations in the opening 
angle. Both changes are shown to have a significant effect on the nematic polymorphism, 
indicating that the lengthening of the molecular arms can help stabilizing the Nb phase. 

The paper is organized as follows. In Sec. II we define model banana-shaped molecules 
and give some computational details. In Sec. HI we introduce relevant order parameters 
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and distribution functions used to identify structures. Results are presented in Sec. IV and 
a short summary is given in Sec. V. 



II. MODEL AND COMPUTATIONAL DETAILS 



We define a banana-shaped Gay-Berne molecule by linking two prolate GB molecules, a 
and b, as shown in Fig. [TJ The positions of the GB centers are given by vectors and r;,. 
In addition, each molecule is parameterized with two internal coordinates: (i) the opening 
angle 7 (O'^ < 7 < 180°) between the unit vectors Ua and u;, entering definition of the GB 
potentials and (ii) the shift c of the GB centers along the molecule's arms. We also introduce 




FIG. 1: Parametrization of model bent-core molecules used in simulations. 

a molecule-fixed orthonormal coordinate system with basis vectors {e'^, /i = 0, 1, 2}, where 
e'2 is taken to belong to the line through endpoints of and r^, e'l is placed along the 
moleular C2-axis, parallel to Uq + u;,, and e'o = e'l x e'2. Details of the parametrization are 
given in Fig. [H 

The intermolecular interaction between a pair {i,j) of the molecules, separated by the 
intermolecular vector Rjj = |[(ra -|- r^-'-') — (ri*^ -I- r^*-*)], is given by: 

f/({et^},{ef^},7.7.,R..) = E E UgbH^ ,i^\r^f), (1) 

a=a,b I3=a,b 

where r[^*^'' = r^"''' — is the intermolecular vector connecting the GB centers (a,/? = 
a, b) of the molecules i and j, respectively, with the vectors u^^"''', uf\ ri?^'' expressed through 



^li^lj^'R-ij and c. Finally, the Ugb term is the Gay-Berne interaction [20|] 
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between the individual arms a, (3 of the molecules i and j, respectively. It reads 

12 



(^0 



6n 



(2) 



where r^*^"''* = ri*^''/''^!*/^ ^'^^ ^^^^ vector and r^*^^'* = Ir^*^"*!. The functions e(ua'', u^"''', r^^"''*) 
and ajuh} , uL^ \ r^^) depend on four GB parameters and p and their explicit form is 

found in j20|]. In addition, each molecule i is allowed to vary its opening angle 7^ about 70 



with a harmonic contribution to the potential energy 

jAi) _ J ^ ^/ l7i - 7o| > Imax , . 

I /i:(7j — 7o) otherwise, 

where /c is the harmonic constant, 70 is the angle for which internal energy has a minimum 
and 7max is the angle of maximal deflection. The case of = 00 corresponds to a rigid 
molecule. 

In our simulations we study the model with the following GB parameters: k = 4, k' = 5.0, 
/i = 1 and z/ = 2. In addition we assume 70 = 140°, 7max = 38°, k = 25,oo and c = 1.7, 
where c is chosen such that a smooth molecular shape is obtained for opening angles in the 
range studied (102° — 178°). Finally, throughout the paper angles are expressed in degrees 



and all other parameters are rendered dimensionless. Their definitions are found e.g. in 13 | 



Memmer [13|] carried out MC simulations for a similar model, but with shorter molecular 
arms. The parameters different from ours were k = 3 and c = 1.0. Only the case with rigid 
molecules (k = 00) was considered. Thus, a comparison of our results with those of Memmer 

n 

|13l ] will allows us to look at the influence of the arms' elongation and of the molecular arms 
fluctuations on the nematic polymorphism. 



All simulations are performed using Metropolis MC simulations in NPT ensemble [21|. 
The system consists of = 1088 molecules placed in a rectangular box with periodic 
boundary conditions and the nearest image convention taken regard of. A sequence of 
configurational states is generated according to the probability density 



= n,^,^(l + 3cos2|) (l -cos^l) cos2|(4-3cos2| 

exp(- ( UtotaiiT, n,^) + PV-NTlnV)/T) / (4) 
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where F, 17, $ denote respectively positions, orientations and opening angles of all molecules, 
V is the volume of the system, Utotai represents the total potential energy and Z'j^pj. denotes 
the partition function. Boltzmann constant, not present in Eq.(jl]), is included into definition 
of the reduced temperature T . The formula (jl]) for the probability density contains square 
root factors coming from the integration over momenta. For rigid molecules the square root 
part gives a constant factor and, hence, cancels out with similar terms in Z'^pj,. 

The MC cycle consisted of probing a new configuration for each molecule followed by 
volume rescaling along a randomly chosen edge of the simulation box. A trial molecular 
configuration was generated by a random move of the molecule's center of mass, a random 
rotation about randomly chosen axis out of the laboratory reference system, a flip of the 
steric dipole, and a random change of the opening angle. In simulations 10% — 50% of new 
molecular configurations and 10% — 60% of system rescalings were accepted. The lower 
acceptance ratios were used for dense, highly ordered states, allowing to partly overcome 
the metastability traps. 

Simulations were carried out under the pressure of P = 3 and with spherical cutoff 
'or the GB potential at the distance r^*^"* = 5.0. The pressure was taken the same as in 
I3I . Simulations were initialized in the isotropic liquid and the temperature was lowered in 
adjustable steps until transition into a highly ordered crystalline smectic phase took place. 
In addition, the crystalline phase found on cooling was heated up all the way back to estimate 
the stability range of the phases found. The equilibration appeared slow and took between 
5 X 10^ and 10^ cycles after which the data were collected for structure analysis once per 10 
cycles over a production run of 10^ cycles. 

Interestingly, the biaxial nematic phase was found on cooling down. To check its stability 
we also melted ground state crystalline structures under higher pressures, from P = 4 to 
P = 10 in step of 1. 



III. ORDER PARAMETERS AND DISTRIBUTION FUNCTIONS 

In order to characterize the long-range orientational order of the molecules the alignment 
tensors were utilized. They are given by 

1 ^ 

Q% = i^{Y.^^lf,l^l-^u)). (5) 

i=l 
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where l^l = e'^^ ■ e^, and {e^j} constitute a base of the laboratory reference system (77, n^u ^ 
{0, 1,2}). Eigenvectors of Q'^^^ correspond to the directors of the phase and the associated 
eigenvalues are the orientational order parameters measured with respect to these directors. 
For isotropic phase all eigenvalues of each tensor Q'^^^ vanish. For uniaxial phase each tensor 
Q"^^^ has one eigenvalue different than the other two, which are equal to each other. In case of 
the biaxial phase all eigenvalues of QJ^^, are, in general, different. However, it might happen 
that one of the tensors be uniaxial [? ]. 

The presence of smectic order was detected by means of the one-particle distribution 
function P^^\ calculated along the main director. More specifically, let the director fi, 
corresponding to the largest eigenvalue of Q^,^, be parallel to 62. Then the smectic ordering 
along n can be examined with the distribution function 

^^^H^) = ^ E - ^ , (6) 

where is the average length of the simulation box along 62, Zi represents coordinate of 
the 2-th molecule along e2, and 5 is the Dirac delta. 

The in-layer structure of the smectic layers was examined with the two-point transverse 
correlation function 

(-^) = ( E E ^(-^ - - i^^u) - 1^^ - , (7) 

where D = 0.5 is the thickness of the slice within which correlations were calculated, V is 
the average volume of the simulation box, and |Rj — Rj|± is the projection of the vector 
linking the origins i, j of the molecule-fixed reference frames on a surface perpendicular to n. 
Owing the small system size we were unable to distinguish between hexatic- and crystalline 
orders and, for this reason, a phase with long-range hexagonal order within layers is referred 
to as a crystal. 

For flexible molecules the one-particle distribution function of opening angles was defined 

as: 



N , 

i=l 



Finally, the relative stability range for a given phase / was examined by calculating the 
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coefficients 



AT^ = (r^-T/)/r/, (9) 

Ap^ = (p^-p{)/p{„ (10) 

where Tj^ {Ph)^ T[ (p{) represent upper and lower limit, respectively, of temperature (den- 
sity), for which the phase / is detected. 



IV. RESULTS AND DISCUSSION 
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Cooling down the system of rigid molecules below T = 4.1 leads to a phase transition from 
isotropic liquid to uniaxial nematic (Figs. [2|3|) . Further lowering of the temperature causes 
stabilization of the biaxial nematic phase, which is found for 2.7 ^ T < 3.0 (Figs. 
Periodic modulation of P^^\z) (Fig. [3]) emerges only below T = 2.8, but for 2.7 ^ T < 
2.8 this modulation is highly irregular, with amplitude not exceeding 0.2. Distinguishable 
smectic order is present only from below T = 2.6 and is accompanied by the hexagonal long- 
range order within smectic layers (see g]_ {r±) shown in Fig. [3]). The phase is antiferroelectric 
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FIG. 2: Eigenvalues Sf^^ of tensors Qf^^ for system of rigid (7 = 140°) and flexible (102° < 7 < 178°) 
molecules obtained on cooling down (V) and heating up (A) the system. The identified phases 
are a) orientationally isotropic, when eigenvalues of are close to zero; b) uniaxial, when each 
tensor Q^j^ has only two different eigenvalues; c) biaxial, when at least two tensors have all 
three eigenvalues different. Note that the biaxial nematic phase appears when the system is cooled 
down. 

with respect to the steric dipoles (Fig. As we cannot distinguish between hexatic and 
crystaUine order we classify this phase as a biaxial, hexagonal and antiferroelectric crystal. 





FIG. 3: Functions g^l\r±) and P^^\z) ob- 
tained for the system of rigid molecules 
(cooling down). For T < 2.7 the system 
gains smectic order with hexagonal long- 
range order within layers. 



FIG. 4: Functions g^l\r_L) and P^^\z) ob- 
tained for the system of flexible molecules 
(on cooling). For T < 2.6 the system gains 
smectic order with hexagonal long-range 
order within layers. 



The crystal melts into uniaxial nematic at T = 3.2. On cooling down the biaxial nematic 
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appeares below T = 3.2; hence without a credible calculation of the free energy it is not 
possible to settle if Nb is really absolutely stable in our simulations (Figs. [2][6l) . Studies 
at higher pressures also do not resolve this issue unambiguously. However, owing that the 
hysteresis affects both measurements it is quite likely that the biaxial nematic might be ab- 
solutely stable in a small temperature interval. Further increase of temperature transforms 
uniaxial nematic into isotropic liquid what takes place at T = 4.4. The coexistence temper- 
ature of isotropic liquid and uniaxial nematic is pinpointed by localizing the discontinuity 
in the density between these phases (Fig. [6]), which yields T = 4.2. 

For the model bent-core molecules with shorter arms of k = 3.0 [l^ no evidence of the 
biaxial nematic phase was reported. Thus we conclude that molecules with more elongated 
arms help to stabilize Nb- In addition, for Memmer's model AT^ = 0.5 and Ap^ = 0.118 
while the same coefficients for our model are AT^ = 0.313 and Ap^ = 0.127. That is the 
elongation of the arms enhances the density range for which nematic phases can be stable. 
Finally, in [isl a helical structure between uniaxial nematic and smectic has been reported. 
Simulations we carried out do not show a formation of any helical order. 




FIG. 5: Left: side view of the crystalline phase obtained for the rigid model at T = 2.6. Steric 
dipoles within neighboring layers have opposite orientations making the crystalline structure "an- 
tiferroelectric" . Right: top view of hexagonal long-range order within smectic layers. 

Cooling down the system of the flexible bent-core molecules leads to the following phase 
transitions: at T = 3.7 the isotropic liquid transforms into the uniaxial nematic and at 
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FIG. 6: Relation between density and 
temperature obtained for the system of 
rigid molecules on cooling down (V) and 
heating up (A). At T = 4.2 isotropic liq- 
uid coexists with uniaxial nematic. Hys- 
teresis between liquid crystalline and crys- 
talline phases (T < 3.2) is also shown. 



FIG. 7: Relation between density and 
temperature obtained for the system of 
flexible molecules on cooling down (V) 
and heating up (A). At T = 3.7 
isotropic liquid coexists with uniaxial ne- 
matic. Shown is also hysteresis between 
liquid crystalline and crystalline phases 
(T < 2.8). 
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FIG. 8: Distribution of opening angles ^^^^^(7) for system of flexible molecules. 

T = 2.7 the nematic phase becomes biaxial (Figs. EH]). For T = 2.6 the one-particle 
distribution function, P'^^^[z\ develops an irregular modulation of magnitude smaller than 
0.2. The ordering becomes periodic along the director for T < 2.5, which is attributed to 
the formation of smectic layers. Within the layers the molecules form again a hexagonal 
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lattice, Fig. HI with snapshots being similar to the one shown in Fig. [51 Hence, the phase 
can again be classified as biaxial, hexagonal and antiferroelectric crystal. 

The crystalline phase of the model melts directly into the uniaxial nematic phase at 
T = 2.8 and, subsequently, the uniaxial nematic transforms into an isotropic liquid at 
T = 3.8. From Fig. [7] we can estimate that the uniaxial nematic and the isotropic liquid 
coexist at T ^ 3.7. Like for the rigid model the presence of hysteresis (Figs. [2][7I) does 
not allow to settle as whether the biaxial nematic is absolutely stable. Note, however, that 
fiuctuations in the opening angle seem to have a destabilizing effect on the biaxial nematic 
phase, as seen in Fig. [21 

For the model with flexible molecules the average opening angle, 7, changes with temper- 
ature (Fig. [8]). At low temperatures, in the crystalline phase, 7 does not differ much from 
7o {i.e. 7 = 140.5° for T = 2.5). As temperature increases the average opening angle first 
increases to 7 = 143° in the uniaxial nematic phase (T = 2.8) to next decrease to 7 = 135° 
in the isotropic hquid phase (T = 3.8). 

The observation of Nb on cooling encouraged us to perform additional simulations at 
higher pressures. In these simulations we limited ourselves to the rigid model of 70 = 140° 
and to integer values of pressure ranging from 4 to 14. In addition, we carried out the 
simulations along the P = 10 isobar for opening angles of 125°, 130° and 135°. In all cases 
we melted the ground state configurations and sought for an interval of stable biaxial nematic 
phase. Apart from the model with the opening angle of 125°, the ground state crystalline 
structures melted directly to the uniaxial nematic phase. For the opening angle of 125° we 
observed a direct melting to the isotropic liquid phase. As a result we cannot confirm that 
the biaxial nematic found on cooling is absolutely stable. 

V. CONCLUSION 

This paper contains results of MC NPT simulations performed for the system of model 
bent-core particles. The arms of V-shaped molecules were built out of two Gay-Berne sites. 
Two variants of the model were studied, with rigid- and with fiuctuating angle between the 
arms. For both models decrease of temperature leads to the sequence of phases: isotropic 
liquid, uniaxial nematic, biaxial nematic and biaxial antiferroelectric hexagonal crystal. On 
heating the ground state configuration the biaxial nematic is superseded by a metastable 
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crystalline phase. However, in contrast to the model with molecules of shorter arms (k = 3.0) 



131], the present systems, especially the rigid one, clearly show a tendency to stabilize Nb- 



In addition, the V-shaped molecules with longer arms extend the range of density in which 



13| was 



the nematic polymorphism can be present. Finally, the helical structure found in 
not detected in the present simulations. 
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